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Abstract 

We study the potential of high-energy linear e + e~ colliders for the production of gluino pairs 
within the Minimal Supersymmetric Standard Model (MSSM). In this model, the process e + e _ — > 
gg is mediated by quark/squark loops, dominantly of the third generation, where the mixing of left- 
and right-handed states can become large. Taking into account realistic beam polarization effects, 
photon and Z°-boson exchange, and current mass exclusion limits, we scan the MSSM parameter 
space for various e + e~ center-of-mass energies to determine the regions, where gluino production 
should be visible. 
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I. INTRODUCTION 



Supersymmetry (SUSY) is generally considered to be one of the most promising exten- 
sions of the Standard Model (SM) of particle physics. Its attractive features include the 
cancellation of quadratic divergences in the Higgs sector, which implies that the soft SUSY 
breaking masses of the (yet unobserved) superpartners of the SM particles can not be much 
greater than the electroweak scale. If SUSY is indeed responsible for the stabilization of 
this scale against the Planck scale, supersymmetric particles should therefore be discovered 
either at Run II of the Fermilab Tevatron |T|, @, §, (§ §] or at the CERN LHC [§, §. In par- 
ticular, the strongly coupling squarks and gluinos should be copiously produced at hadron 
colliders and lead to first measurements of their masses and production cross sections |§. 
Precision measurements of masses, mixings, quantum numbers, and couplings must, how- 
ever, be performed in the clean environment of a future linear e + e~ collider because of the 
large hadronic SM background and theoretical scale and parton density uncertainties at the 
Fermilab Tevatron and CERN LHC. For example, in e + e~ annihilation the center-of-mass 
energy of the collision is exactly known, and threshold energy scans allow for a precise mass 
determination of pair-produced SUSY particles. It will then be possible to establish whether 
the masses and couplings of the electroweak gauginos and of the gluino are indeed related, as 
expected. A global analysis should ultimately lead to a reconstruction of the SUSY breaking 
model and its parameters. Along these lines, detailed studies have recently been performed 
for squarks, sleptons, charginos, and neutralinos ||, but not for gluinos, the reason being 
that gluino pairs are produced at e + e~ colliders only at the one-loop level, while all other 
sparticles are produced at tree level. At tree level, gluinos can be produced in pairs only 
in association with two quarks fT(| , or they are produced singly in association with a quark 
and a squark JTI], [T2|]. Both processes result in multi-jet final states, where phase space is 
limited and gluinos may be hard to isolate. 

In the Minimal Supersymmetric Standard Model (MSSM) |L3[ [14}] , the exclusive produc- 
tion of gluino pairs in e + e~ annihilation is mediated by s-channel photons and Z°-bosons, 
which couple to the gluinos via triangular quark and squark loops. In earlier studies of this 
process only the very low center-of-mass energy region (y^ = 20 GeV) with pure photon 
exchange and no squark mixing [[DJ or Z°-boson decays into light (nig < raz/ty JX6|, [T7| 
and very light (rrig = 3 ... 5 GeV) [18|, [19| gluinos have been considered. Some authors have 



2 



presented results only for nig = GeV and outdated top quark masses of 20... 50 GeV 
|T6|, |I9||, while others have neglected the mixing of left- and right-handed squark interaction 



eigenstates into light and heavy mass eigenstates [|I^, [L6|, pH , which turns out to control the 
production cross section to a large extent. 

It is the aim of this Article to study the potential of high-energy linear e + e~ colliders 
for the production of gluino pairs within the MSSM. Taking into account realistic beam 
polarization effects, photon and Z°-boson exchange, and current mass exclusion limits, we 
scan the MSSM parameter space for various e + e~ center-of-mass energies to determine the 
regions, where gluino production should be visible. Furthermore we clarify the theoretical 
questions of the relative sign between the two contributing triangular Feynman diagrams, 
of the possible presence of an axial vector anomaly, and the conditions for vanishing cross 
sections - three related issues, which have so far been under debate in the literature. The 
remainder of this Paper is organized as follows: In Sec. II, we present our analytical results 
and compare them with existing results in the literature. Various numerical cross-sections 
for gluino pair production at future high-energy linear e + e~ colliders are computed and 
discussed in Sec. [TTT], and Sec. [TV] contains our Conclusions. Our conventions for squark 
mixing are defined in App. and a summary of all relevant Feynman rules is given in 
App. 0. 



II. ANALYTICAL RESULTS 

The scattering process 

e~(pi, Ai)e + (>2, A 2 ) -> g(h)g(k 2 ) (1) 

with incoming electron/positron momenta px % and helicities A12 and outgoing gluino mo- 
menta k\p proceeds through the two Feynman diagrams A and B in Fig. [I] with s-channel 
photon and Z°-boson exchange and triangular quark and squark loops. Higgs boson ex- 
change is not considered due to the negligibly small electron Yukawa coupling, but it could 
well be relevant at muon colliders. The process occurs only at the one-loop level, since the 
gluino as the superpartner of the gauge boson of the strong interaction couples neither di- 
rectly to leptons nor to electroweak gauge bosons. Taking into account chiral squark mixing 
(see App. [A]) and using the Feynman rules in App. [B], we decompose the corresponding 
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FIG. 1: Feynman diagrams for gluino pair production in electron-positron annihilation. The 
exchanged photons and Z°-bosons couple to the produced gluinos through triangular qqqi (A) and 
QiQjQ (B) loops with flavor flow in both directions. 



scattering amplitude 



into the lepton current 
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the photon and Z°-boson propagators 
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which depend on the squared center-of-mass energy s — (pi + p 2 ) 2 , and the gluino current 
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q L i 
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and 



m q ){2q - k x + k 2 ) l 



(2tt) d (q 2 — m 2 + irj) \{q — k\) 2 — mj + irj] [(q + k 2 ) 2 — mf + irj] 



(7) 
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and the quark flavor q flowing both ways in the corresponding diagrams A and B, so that 
= A % V{v v q -> -v\), Bf Y = -B l J> v , and V = C^C' 1 = T |20|. 



Eq. (El) can be simplified with the Dirac equation and the anti-commutation relations for 
Dirac matrices, and the tensor loop integrals in Eqs. (0) and (|7|) can be expressed through 
the standard coefficient functions C^n of the metric tensor g^ u and tensors constructed from 
the outgoing gluino momenta ki tlx and k^^ [^U], PT] . The gluino current then reduces to 



Gl = ie-^ —u(k 2 )w 5 v(h) £(< + B v q ) (8) 



27T 2 



with 



A q=J2 [ C o\ m l a qiV - m l a qiv + ^m q m q a qW ) + CfAm q (m q a qiV - m~ g a+ v ) (9) 

i 

+Cg(2 - -D) a+ v - Cpmla+y + Cf 2 (s - 2mf) a+ v , 
B V q = E^2b qijV: (10) 



where Cg ;) = Cjfc(j)( m l> s, mj, m|, m 2 q , m 2 q ) and = C 00 (m 2 g} s, m|, mj, m|, m|) are mas- 
sive (infrared-finite) three-point functions and Cf = Cf* and C'} = CI2 in diagram A. 

n ± _ ,y (Qi qi* _ qi qi*\ _l v 

a qiV — U q K^il^il t ->i2'->i2 ) ^ a q > 



a<7iv = (^(SfiSg + Sf 2 Sf*) , and (11) 

UqijV — D il D jl L q D i2°j2 L q 



are combinations of vector (v Y ), axial vector (a Y ), and derivative couplings (T q ^ ), and 
elements of the squark mixing matrix S. Pairs of identical (Majorana) gluinos are therefore 
produced by a parity violating axial vector coupling induced by mass differences between the 
chiral squarks and the axial vector coupling a q of the Z°-boson. The (mass-independent) ul- 
traviolet singularities contained in the Coo-functions cancel among A Y and B Y in D = 4 — 2e 
dimensions. As we have checked explicitly (even for complex squark mixing matrices), adding 
the two amplitudes induces not only a cancellation of the ultraviolet singularities and of the 
logarithmic dependence on the scale parameter /z introduced in Eqs. (H) and (|7|) in order 
to preserve the mass dimension of the loop integrals, but also a destructive interference of 
the finite remainders. This happens separately for each weak isospin partner, as is to be 
expected for triangular loop diagrams involving one axial vector and two scalar (not vector) 
couplings and no closed fermion loop. 



The (finite) total cross section for incoming electrons/positrons with helicities Ai,2 
is then 



with 
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= (« 2 + <a™){l - 4A X A 2 ) - (^ar 2 + v^){2X x - 2A 2 



±1/2 



;i2) 



(13) 



color factor iVc = 3, and gluino velocity /3 — yl — 4m|/s, which contains the expected fac- 
tors of /3 3 and s for P-wave production of two spin-1/2 Majorana fermions. The distribution 
in the center-of-mass scattering angle 9, 



( S ) = ^(l + cos 2 ^ AlA2 ( S ), 



(14) 



is independent of the gluino mass and has to be integrated over just one hemisphere, since 



the two final state particles are identical JL8J, p22| . As a consequence, the forward-backward 
asymmetry vanishes for Majorana fermions, but not for Dirac fermions. 
Our result for diagram A agrees with the unpolarized result 

= 7 £ (J A 1 A 2 (s) 
4 Ai |2 =±l/2 

in Eq. (4.5) of Kileng and Osland, if we identify |18, ^ 



(15) 



ni' 1 
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°00 
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_P00 
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(16) 



and reverse the sign of b q to account for opposite conventions of squark mass eigenstates. 
However, our result for diagram B disagrees in sign with Eq. (4.5) of Kileng and Osland, if 
we identify |1| ^3 



c: 



00 



-Gijg/2. 



(17) 



If the sign of diagram B is reversed, the ultraviolet singularities cancel only after adding the 
contributions from the two weak isospin partners with opposite values of T 3 . We trace this 



sign discrepancy to the Feynman rules employed in Ref . ]18| , which exhibit a relative minus 
sign to those in Ref. |14[] for the Z°-boson coupling to quarks, but not to squarks, whereas 
our Feynman rules (see App. |B|) agree with those in Ref. []I4J in the limit of no squark mixing. 
Except for the relative sign of diagrams A and B and in the limit of vanishing gluino mass, 



we also find agreement with Djouadi and Drees [B5|. We confirm, however, the relative sign 
for diagrams A and B of Campbell, Scott, and Sundaresan jT7|, who found (for non-mixing 
chiral squarks of different mass) that the ultraviolet singularities cancel separately for each 
weak isospin partner, and that there is no anomaly. This had also been claimed previously 



by Kane and Rolnick |16[ for chiral squarks of equal mass. In their limit, the cross section 



depends only on the weak isospin and not on the charge of the (s) quarks [16], and the 
contribution of the photon vanishes [15] . For the contribution of the Z°-boson to vanish, we 



must have 18 



1. mass degeneracy in each quark isospin doublet, = m u etc., 

2. mass degeneracy in each squark isospin doublet, m? = mj 2 = m fil = etc., 



which contradicts the condition m q = rriq found by Kane and Rolnick [IB]. Condition (1) is 
violated most strongly for the third generation, as is condition (2) for most SUSY breaking 
models. 



III. NUMERICAL RESULTS 

The analytical results presented in the previous Section have been obtained in two inde- 
pendent analytical calculations. They have been implemented in compact Fortran computer 
codes, which depend on the LoopTools/FF library ^] for the evaluation of the mas- 
sive tensor three-point functions. As a third independent cross-check, we have recalculated 
the production of gluino pairs in e + e~ annihilation with the computer algebra program 
FeynArts/FormCalc [^!| and found numerical agreement up to 15 digits. 

Our calculations involve various masses and couplings of SM particles, for which we 



use the most up-to-date values from the 2002 Review of the Particle Data Group |27| . In 
particular, we evaluate the electromagnetic fine structure constant a(mz) = 1/127.934 at 
the mass of the -Z^-boson, mz = 91.1876 GeV, and calculate the weak mixing angle 8w 



from the tree-level expression sin 2 6 W = 1 — m^/ml with mw = 80.423 GeV. Among the 



fermion masses, only the one of the top quark, m t = 174.3 GeV, plays a significant role due 
to its large splitting from the bottom quark mass, m& = 4.7 GeV, while the latter and the 
charm quark mass, m c = 1.5 GeV, could have been neglected like those of the three light 
quarks and of the electron/positron. The strong coupling constant is evaluated at the gluino 
mass scale from the one-loop expression with five active flavors and Alo~ 5 = 83.76 MeV, 
corresponding to a s (mz) = 0.1172. A variation of the renormalization scale by a factor of 
four about the gluino mass results in a cross section uncertainty of about ±25 %. Like the 
heavy top quark, all SUSY particles have been decoupled from the running of the strong 
coupling constant. 

We work in the framework of the MSSM with conserved R- (matter-) parity, which rep- 
resents the simplest phenomenologically viable model, but which is still sufficiently general 
to not depend on a specific SUSY breaking mechanism. Models with broken i?-parity are 
severely restricted by the non-observation of proton decay, which would violate both baryon 
and lepton number conservation. We do not consider light gluino mass windows, on which 
the literature has focused so far and which may or may not be excluded from searches at fixed 



target and collider experiments [2~7|. Instead, we adopt the current mass limit m~ g > 200 GeV 



from the CDF []28| and DO [^9| searches in the jets with missing energy channel, relevant 
for non-mixing squark masses of rriq > 325 GeV and tan (5 = 3. Values for the ratio of the 
Higgs vacuum expectation values, tan (3, below 2.4 are already excluded by the CERN LEP 
experiments |30|| . If not stated otherwise, we will present unpolarized cross sections for a 
t/s = 500 GeV linear e + e~ collider like DESY TESLA, gluino masses of rrig = 200 GeV, and 
squark masses rriq ~ rriQ — TO5 = rriy = m§ = = = rrif = msusY = 325 GeV. We 
will consider two cases of large squark mass splittings: I.) On the one hand, the masses of 
the superpartners of left- and right-handed quarks need not be equal to each other. In this 
scenario we will vary the right-handed up-type squark mass parameters rrifj q f between 200 
and 1500 GeV. II.) On the other hand, the superpartners of the heavy quarks can mix into 
light and heavy mass eigenstates (see App. |A|). This alternative is restricted by the CERN 
LEP limits on the light top and bottom squark masses, > 100 GeV and m~ bi > 99 GeV 
PHI , and on SUSY one-loop contributions [^, |33|, [54J to the p-parameter, Psusy < 0.0012 
| 27|| . In this case we assume the maximally allowed top squark mixing with 9$ = 45.2°, 
= 110 GeV, and mi 2 = 506 GeV, which can be generated by choosing appropriate 
values for the Higgs mass parameter, // = —500 GeV, and the trilinear top squark coupling, 
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At = 534 GeV. For small values of tan /3, mixing in the bottom squark sector remains small, 
and we take 6~ b = 0°. Although the absolute magnitude of the cross section depends strongly 
on the gluino mass and collider energy, the relative importance of the different contributions 
is very similar also for higher gluino masses and collider energies. 

First we examine the conditions found in Section II for vanishing of the photon and 
.Z^-boson contributions, restricting ourselves to the third generation. Since we expect the 
photon contribution to cancel for equal left- and right-handed squark masses, we vary the 
right-handed top squark mass parameter, mf ~ m t R -> between 200 and 1500 GeV, but 
keep mi L ~ nig = mp, = mgusY = 325 GeV fixed (case I), since top and bottom squarks 
generally interfere destructively due to their opposite charge and weak isospin quantum 
numbers. As can be seen from Fig. |2], the photon contribution cancels indeed for rrif ~ 
mi R = mi L ~ msusY- This is due to the fact that for photons a qil = and 6 g i2 7 = 
b q 2i~/ = in Eq. (|TTD, while unitarity of the squark mixing matrix leads to a^ l7 = — a^ 7 
and 6 9 n 7 = — & 9 22 7 - Therefore, the photon contributions cancels for all flavors q with equal 
squark masses. Due to their charge, top (s)quarks contribute four times as much as bottom 
(s) quarks, whose contribution is even more suppressed by the condition ~ m~ b . The 
Z°-boson contribution can never cancel, since m t ^> mj, and therefore it depends only 
weakly on mf, but it can become minimal for mf ~ m t R = m t L = m l R = m l L — m susY- 
As mf gets significantly larger (or smaller) than mgusY; the photon contribution starts to 
dominate over the Z°-boson contribution. 

If only mf differs from msusY, the third generation contributes almost 100% to the total 
cross section. However, if m^ = m^ = mf are varied simultaneously, all three generations 
contribute to the total cross section, which can therefore become significantly larger. This 
is shown in Fig. |3], where (s) quark loop contributions from all three generations have been 
taken into account. 

When mjj = m^ = mf = msusY and large mass splittings are generated only by mixing in 
the top squark sector (case II), photon contributions are suppressed by more than two orders 
of magnitude. Fig. ^ shows that the Z°-boson contributions from top and bottom squarks 
interfere destructively due to opposite values of their weak isospin quantum numbers, except 
for 9i ~ 45.2°, where the imaginary parts of the amplitudes interfere constructively. It is 
therefore advantageous to keep the bottom squark mass splitting small. As is also evident 
from Fig. ^, mixing in the bottom squark sector is of little importance. Note that the central 
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m~ Dependence of y, Z Contributions 
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FIG. 2: Dependence of the photon and Z°-boson contributions to the process e + e~ — > gg on the 
right-handed top squark mass parameter rrif. The photon contribution (dashed curve) is dominated 
by top (s)quarks and cancels for = rn^ R . The Z°-boson contribution from top (dotted curve) 
and bottom squarks (dot-dashed curve) interferes constructively with the photon contribution (full 
curve) . 

region with maximal top/bottom squark mixing is excluded by the CERN LEP limits on 
mi x , mg , and the p-parameter. 

When mgusY and the diagonal elements of the squark mixing matrix (see App. [AJ) become 
much larger than the quark masses and the off-diagonal elements of the matrix, the role of 
squark mixing is expected to be reduced. This is confirmed numerically in Fig. |5|, where 
the dependence of the gluino production cross section on msusY is shown for the cases of 
maximal and vanishing top squark mixing. Squarks from the first two generations contribute 

10 




0.25 
0.225 
0.2 
0.175 
0.15 
^0.125 
0.1 
0.075 
0.05 
0.025 



Mass Dependence of q Loop Contributions 
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FIG. 3: Dependence of the q loop contributions to the process e + e~ —* gg on the right-handed up- 
type squark mass parameter m^y = m^, = rrif. When these mass parameters differ simultaneously 
from msusY = 325 GeV (full curve) or msusY = 1000 GeV (dashed curve), all three (s)quark 
generations contribute significantly to the total cross section, so that it becomes much larger than 
in the case where only rrif is varied (dotted curve). 

at most 10% at low msusY and are otherwise strongly suppressed. 

At future linear e + e~ colliders it will be possible to obtain relatively high degrees of 
polarization, i.e. about 80% for electrons and 60% for positrons 0. In Fig. |6| we there- 
fore investigate the effect of choosing different electron/positron polarizations on the gluino 

pair production process, including contributions from all (s) quarks. Since the ++ and 

helicity amplitudes vanish for both photons and Z°-bosons, we only show the squares of 
the remaining H — and — h amplitudes, which coincide for photons, but not for Z°-bosons. 
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Mixing Angle Dependence of t, b Loop Contributions 
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FIG. 4: Mixing angle dependence of the t (dashed) and b (dotted) loop contributions to the process 
e + e~ — > gg, which interfere destructively (full curve), except for Of ~ 45.2°, where the imaginary 
parts of the amplitudes interfere constructively. Mixing in the b sector (dot-dashed curve) enhances 
the cross section only slightly. 

The unpolarized cross section falls short of the polarized ones, so that a high degree of 
polarization is clearly desirable. 

With the realistic degrees of polarization mentioned above, we show in Fig. [7| a scan in 
the center-of-mass energy of a future e + e~ collider for various gluino masses and maximal 
top squark mixing (case II). The cross section rises rather slowly due to the factor /3 3 in 
Eq. (|i~2D for P-wave production of the gluino pairs. For m- g = 200 GeV we observe an 
interesting second maximum, which arises from the intermediate squark pair resonance at 
^Js = 2 msusY = 650 GeV. At threshold, the cross section depends strongly on the gluino 
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Mass Dependence of q Loop Contributions 
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FIG. 5: Squark mass dependence of the loop contributions from third generation squarks with 
maximal (full curve) and vanishing mixing (dashed curve). The contributions from the first two 
generations (dotted curve) are highly suppressed. 

mass and is largest for rrig = 200 GeV, which we consider to be the lowest experimentally 
allowed value. It drops fast with trig, so that for rrtg > 500 GeV no events at colliders with 
luminosities of 1000 fb _1 per year can be expected, irrespective of their energy. Smaller 
squark mixing (cf. Fig. ^) or larger values of ?«susy (cf. Fig. [5]) will reduce the cross section 
even further. Far above threshold, it drops off like 1/s and becomes independent of the 
gluino mass. 

The slow rise of the cross section can be observed even better in Fig. || where the 
sensitivity of a ^Js = 500 GeV collider like DESY TESLA to gluino masses around 200 GeV 
has been plotted. For the CERN LHC experiments, a precision of ±30 ... 60 (12 ... 25) GeV 
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Polarisation Effects in y, Z Contributions 



-l 




Vs~[GeV] 

FIG. 6: Center-of-mass energy dependence of the process e + e~ — > gg for unpolarized (full curve) 
and polarized (dashed and dotted curves) incoming electrons/positrons and maximal top squark 
mixing. The photon contribution (dot-dashed curve) is suppressed by more than two orders of 
magnitude. 

is expected for gluino masses of 540 (1004) GeV H, [7J. If the masses and mixing angle(s) 
of the top (and bottom) squarks are known, a precision of ±5... 10 GeV can be achieved 
at DESY TESLA for m g = 200 GeV and maximal top squark mixing with an integrated 
luminosity of 100 fb _1 per center-of-mass energy point. 

A center-of-mass energy scan for the scenario with no squark mixing, but large left- 
/right-handed squark splitting (case I) is shown in Fig. |9| for light (full and dashed curves) 
and heavy (dotted and dot-dashed curves) gluino masses. Since the photon contributes 
now significantly to the cross section, it proves to be advantageous to choose the lepton 
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e (80%-) e + (60%+) -> g g 
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FIG. 7: Center-of-mass energy dependence of the polarized e + e 
gluino masses and maximal top squark mixing. 



gg cross section for various 



polarization such that the Z°-boson interferes constructively with the photon, even though 
it is by itself slightly smaller than for the opposite choice. Since all three generations add 
now to the cross section, it can become almost an order of magnitude larger than in the 
mixing scenario (case II), and even gluino masses of 1 TeV may be observable at a multi- 
TeV collider like CERN CLIC. However, also here the cross section drops sharply when the 
squark mass splitting is reduced from 1500 TeV (full and dotted curves) to values close to 
^susy (dashed and dot-dashed curves). 
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FIG. 8: Sensitivity of the polarized e + e~ — * gg cross section to the gluino mass m g for maximal top 
squark mixing. The central values and statistical error bars of the data points have been calculated 
assuming m g = 200 GeV and a luminosity of 100 fb _1 per center-of-mass energy point. 

IV. CONCLUSIONS 

The purpose of this Paper has been two- fold: First, we have resolved a long-standing dis- 
crepancy in the literature about the relative sign of the quark and squark loop contributions 
to the production of gluino pairs in e + e~ annihilation. We confirm the result of two older 
papers that the divergence cancels for each squark flavor separately and not between weak 
isospin partners [|T6| , [TTj] and trace the sign problem in one case to the Feynman rules em- 
ployed in the corresponding calculation fll8f . Our results rely on two completely independent 
analytical calculations and one computer algebra calculation. 
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FIG. 9: Center-of-mass energy dependence of the polarized e + e — > gg cross section for various 
gluino masses and mass splittings between left- and right-handed up-type squarks. 

Second, we have investigated the prospects for precision measurements of gluino proper- 
ties, such as its mass or its Majorana fermion nature, at future linear e + e" colliders. We 
have taken into account realistic beam polarization effects, photon and Z°-boson exchange, 
and current mass exclusion limits. Previously, only light gluinos at center-of-mass ener- 
gies up to the Z°-boson mass had been investigated. Within the general framework of the 
MSSM, we have concentrated on two scenarios of large left-/right-handed up-type squark 
mass splitting and large top squark mixing, which produce promisingly large cross sections 
for gluino masses up to 500 GeV or even 1 TeV. Gluino masses of 200 GeV can then be 
measured with a precision of about 5 GeV in center-of-mass energy scans with luminosities 
of 100 fb _1 /point. However, when both the left-/right-handed squark mass splitting and 
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the squark mixing remain small, gluino pair production in e + e annihilation will be hard to 
observe, even with luminosities of 1000 fb~ 1 /year. 
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APPENDIX A: SQUARK MIXING 

The (generally complex) soft SUSY-breaking terms A q of the trilinear Higgs-squark- 
squark interaction and the (also generally complex) off-diagonal Higgs mass parameter /i 
in the MSSM Lagrangian induce mixings of the left- and right-handed squark eigenstates 



qi,,R of the electroweak interaction into mass eigenstates q\^- The squark mass matrix [14], |35 

M 2 



' m 2 LL + m\ m q m* LR 



\ m q m LR m\ R + m\ 



(Al) 



with 



m 



\ L = (T 9 3 -e q sm 2 9 w )m 2 z cos2p + ml, (A2) 



„ 22 mf, for up — type squarks, 

m RR = eg sin Qw m z cos 2/5 + \ (A3) 

I m 2 ^ for down — type squarks, 

{cot (3 for up — type squarks 
(A4) 
tan (3 for down — type squarks 

is diagonalized by a unitary matrix S, SAi 2 S Ji = diag(mf,m|), and has the squared mass 
eigenvalues 

m h = m l + \ ( m LL + m RR T V ( m LL - m R R ) 2 + 4m 2 q \m LR \ 2 ) . (A5) 
For real values of ttilr, the squark mixing angle 9 q , < 9 q < n/2, in 

5=| «"*• "^O with f*|=S| " | (A6) 
- sin 6 q cos 6 q J \q 2 
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can be obtained from 



tan 20,= 2 ™« mL « . (A7) 



m LL ~ m RR 



If itilr is complex, one may first choose a suitable phase rotation q' R = e^g^ to make the 
mass matrix real and then diagonalize it for q^ and q' R . tan (3 = v u j Vd is the (real) ratio of the 
vacuum expectation values of the two Higgs fields, which couple to the up- and down-type 
(s)quarks. The weak isospin quantum numbers for left-handed up- and down-type (s)quarks 
with hypercharge Y q = 1/3 are Tjj = {+1/2, -1/2}, whereas Y q = {4/3, -2/3} and Tjj = 
for right-handed (s)quarks, and their fractional electromagnetic charges are e q = T g 3 + Y q /2. 
The soft SUSY-breaking mass terms for left- and right-handed squarks are rrig and m^, m^, 
respectively, and rriz is the mass of the neutral electroweak gauge boson Z°. 

APPENDIX B: FEYNMAN RULES 

Denoting squark mass eigenstates by i, j, Lorentz indices by /i, z/, and color indices 
of the fundamental (adjoint) representation of the color symmetry group SU(3) by /, m, ... 
(a, b, ...), we obtain the following propagators in Feynman-gauge: 

V = Y ' Z ° -ia^ 
H wvwyww* v ]9 m 2 = {Q) m || (B1) 

p p 2 — rriy + it] 

l . -2— . m Kt> + rn q )6i m 



j,m 



p2 _ m 2 _|_ ^ 



idind, 



(B2) 



(B3) 



p p z — m~ + il] 

Dirac fermions carry an arrow, which indicates the fermion number flow, whereas Majorana 
fermions, such as gluinos, do not carry arrows. An additional arrow is depicted next to all 
fermion lines in order to obtain a unique orientation of the fermion flow, which is evaluated 
according to the rules in f36 |. The interaction vertices are given by [14| 




\v£yr v v \ v J = e f a / = ° 

iej,(v f - a fl5 )S lm { ^ T! _ 2efSin28w ^ TJ (B4) 

/ 2 sin 9w cos 9w f 2 sin 9w cos 9w 
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n -*e(p+pVfXn (B5) 

\\ P Y ij ' Z ° = (r 9 3 -e g sin 2 e w )S jl S* 1 -e q sin 2 ^S j2 S' 2 

q S \ I. 1 sin 6*14/ cos 



m 

iTl, = -iy/2g s T? m {S iX P L - S l2 P R ) (B6) 

i.i 

m 

iT% 2 = -iyfig.Thi&PR ~ S* 2 P L ) (B7) 
X 

9 \ 

s y 

The Feynman rules in Eq. ( |B4]) apply to photon and Z°-boson interactions with quarks 
and charged leptons. The latter carry electromagnetic charge = — 1 and weak isospin 
Tf = —1/2 (left-handed) and (right-handed), but no color (5i m — > 1). The gauge couplings 
g and g' of the weak isospin and hypercharge symmetries SU(2)l and U(l)y have been 
expressed in terms of the electromagnetic coupling e = g sin 9w and the sine and cosine of 
the weak mixing angle sin Ow/ cos 9w = g'/g- The gluino-quark-squark vertices depend on 
the generators of the SU(3) color symmetry group, Tfi^, and on the Yukawa coupling g s , 
which is identical to the strong gauge coupling g s in leading order, and on the squark mixing 
matrix S, but not on the orientation of the fermion flow. 
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